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We calculate the primordial black hole (PBH) mass spectrum produced from a collapse of the 
primordial density fluctuations in the early Universe using, as an input, several theoretical models 
giving the curvature perturbation power spectra Vn(k) with large (~ 10 -2 -r 10 -1 ) values at some 
scale of comoving wave numbers k. In the calculation we take into account the explicit dependence 
of gravitational (Bardeen) potential on time. Using the PBH mass spectra, we further calculate the 
neutrino and photon energy spectra in extragalactic space from evaporation of light PBHs, and the 
energy density fraction contained in PBHs today (for heavier PBHs). We obtain the constraints on 
the model parameters using available experimental data (including data on neutrino and photon 
cosmic backgrounds). We briefly discuss the possibility that the observed 511 keV line from the 
Galactic center is produced by annihilation of positrons evaporated by PBHs. 

PACS numbers: 98.80.-k, 04.70.-s 
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I. INTRODUCTION 

It is well known that for a sufficient production of 
primordial black holes (PBHs) in the early Universe 
IE & H, 13 the spectrum of the density perturbations 
set down by inflation must be "blue", i.e., it must have 
more power on small scales. This implies that the spec- 
tral index of the scalar perturbations must be larger than 
L in strong contradiction with the latest WMAP results 
[H> B • In particular, standard inflationary models of 
hybrid type, in which the inflaton is trapped in a local 
minimum of the potential and which predict blue spectra 
seem to be excluded as a possible source of PBHs. 

Such a conclusion is correct, however, only in rather 
special case: namely, it is based on the prediction of slow- 
roll single field inflationary scenario, according to which 
the power spectrum of curvature perturbations is nearly 
scale-invariant, i.e., the spectral index n is close to unity 
and the variation in the spectral index dn/dlog k is small. 

Although a prediction of the approximate scale in- 
variance of the primordial power spectrum is a neces- 
sary requirement to any inflationary model, some de- 
viations from pure scale invariance are consistent with 
the observational data. These deviations are described 
by adding localized features to the primordial spectrum 
(see, e.g., [|| and references therein) and/or by intro- 
ducing spectral features modifying a single power law. 
Models with such peculiarities (sometimes called broken- 
scale-invariant (BSI) models) were proposed, in main as- 
pects, in eighties [|, M, M, E Gl M- Such models 
generally include, in addition to the usual inflaton field, 
other scalar fields driving successive stages of inflation 
and triggering phase transitions. 

Evidently, the BSI models of inflation could predict, 
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generically, the essential production of primordial black 
holes at small and medium scales. In particular, in [l5| 
the inflationary model with a steplike power spectrum 
based on a nonanalyticity in the inflaton potential was 
proposed and in subsequent works the inflationary po- 
tentials of such kind were used for making predictions of 
PBH production [IE HOI]- 

A step feature in the potential is an effective field the- 
ory description of a phase transition, so, in a more re- 
alistic approach an inflationary model should certainly 
involve more than one scalar field, and such rapid phe- 
nomena as phase transitions [19| . Second order phase 
transitions during inflationary expansion had been first 
considered in [TTI. Il3l] in models with two scalar fields. In 
scenarios of such type, during a short stage, correspond- 
ing to the beginning of a phase transition, the mass of the 
trigger field becomes negative, and adiabatic perturba- 
tions are exponentially amplified resulting in the forma- 
tion of a narrow spike in the primordial spectrum, and, as 
a consequence, in a copious production of PBHs [20l l2l| . 
There are many multiple field scenarios predicting the 
existence of spike- or bumplike features in the primor- 
dial spectrum (e.g., supersymmetric double hybrid mod- 
els [22j . multiple inflation models based on supergravity 
[23| , etc) . Some of these models are specially constructed 
to predict efficient PBH production [24|, |25| . 

Further, large curvature perturbations (especially on 
subhorizon scales) leading to features in the primordial 
spectrum and possible PBH production, can arise in mul- 
tiple field scenarios at the end of inflation (during pre- 
heating era) or between two consecutive stages of infla- 
tion, as a result of parametric resonance or tachyonic in- 
stability. Estimates of PBH production possibilities had 
been obtained in [5E [U [H[ for two-field models and in 
[29l ] for models with a self-interacting scalar field. The 
more complicated model (contai ning three fields), based 
on supergravity, is considered in [2a ]. 

An existence of the narrow spikes in the primordial 
spectrum is possible not only in multiple field inflation- 
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ary scenarios. Such a feature can, in principle, exist even 
in single field models (see, e.g., [30,[3f(). If, in particular, 
the inflationary potential has an unstable maximum at 
origin (e.g., the double- well potential) then, with some 
fine-tuning of parameters and initial conditions, the in- 
flation process may have two stages, with a temporary 
stay at the maximum, that may lead to the correspond- 
ing peak in the primordial spectrum and, depending on 
the amplitude of the peak, to the PBH production. 

Another example of an inflationary model predicting 
large amplitudes of the density perturbations at small 
scales is the inflationary model with the running mass 
potential. More generally, large amplitudes of the pri- 
mordial curvature perturbation spectrum at small scales 
are possible in models of "hilltop hybrid inflation" [32| . 
Potentials of these models have concave-downward form 
at cosmological scales (corresponding to red primordial 
spectrum as required by observations), but can be much 
flatter near the end of inflation. Such a modification 
of a standard hybrid inflation scenario is discussed, e.g., 
in [33|, [34j . The running mass potential belongs to the 
class of hilltop potentials; in models with such potentials, 
due to spectral index scale-dependence, the amplitude of 
the perturbation spectrum at small scales depends on 
the value of dn/d\ogk at cosmological scales 32[. It 
is shown in our previous work [3l| that possibilities of 
noticeable production of PBHs in running mass model 
are still open in spite of the rather severe experimental 
constraints [35|, [3(| on the value of the spectral index 
running. 

In this paper we limit ourselves to the situations where 
PBHs form from the density perturbations, induced by 
quantum vacuum fluctuations during inflationary expan- 
sion. The details of the PBH formation had been studied 
in [53, H^|, the astrophysical and cosmological constraints 
on the PBH density had been obtained in many subse- 
quent works (see, e.g., the recent review [^). The or- 
der of magnitude of the corresponding constraint on the 
value of the density perturbation amplitude (for the PBH 
mass region 10 11 — I0 18 g, which we are interested in) is 
well known [39], but, if the primordial spectrum contains 
the peaklike feature, the concrete value of the PBH con- 
straint clearly depends on the parameters characterizing 
the form of this feature (in particular, on the width of 
the peak). Such an information may be rather useful for 
the model makers. 

One must note that the formation of PBHs in models 
with the primordial spectrum having the features had 
been studied earlier, in works [iiJIHj]. In the case of 
the spectrum that is sharply peaked on a single mass scale 
was considered. The /3(M) function (the probability of 
a region of mass M to form a PBH) was calculated and 
it was shown that the form of this function depends on 
the type of the gravitational collapse. In :4f| two cases 
were studied: a pure step in the primordial spectrum 
and the spectrum produced in an inflationary model [l5T ] 
with a jump in a first derivative of the inflaton potential 
at some scale. It was shown that the corresponding (5- 



functions have the pronounced bumps and, in connection 
with these bumps, authors of [4I[ discuss the possibility 
that PBHs can be a significant part of dark matter. 

In the present work we reconsider the problem of con- 
straining the power spectrum of the primordial fluctu- 
ations (with accent on the peaklike features) calculat- 
ing the process of the formation of PBHs having small 
masses (~ 10 11 — 10 18 g). Products of evaporation of these 
PBHs contribute, in particular, to extragalactic photon 
and neutrino diffuse backgrounds (which are measured 
experimentally) , and the constraints are calculated using 
the standard procedure [42| . We do not calculate the 
/3-functions preferring to constrain the power spectrum 
directly. The concrete constraints are obtained for two 
special cases (the power spectrum with a peak and the 
spectrum of the running mass model), although the gen- 
eral formulas of Sec. [TT1 and HTT1 can be used for any form 
of the spectrum. 

At the end of the paper we use our approach for 
a checking, once more, the idea proposed in [H, Hij . 
namely, we study the possibility that evaporating PBHs 
(having mass spectrum calculated with formulas of Sec. 
HI)) cluster in the Galactic center and produce the ob- 
served 511 keV photon line. 

The plan of the paper is as follows. In Sec. |TT] we 
present the general formalism of the PBH mass spectrum 
calculation with taking into account the explicit time de- 
pendence of the gravitational potential. The results of 
such a calculation for the concrete example, in which the 
curvature perturbation spectrum has a sharp peak, are 
given. In Sec. Mil we give some details of the calcula- 
tion of the neutrino and the photon extragalactic diffuse 
spectra from the PBH evaporations. In Sec. IIVI and Sec. 
IVl the constraints on the curvature spectrum from PBHs 
are given, for the model of V-jz(k) with a peak and for 
the running mass model, respectively. In Sec. IVII we 
consider the possibility of the explaining the observed 
511 keV photon line from the Galactic center by PBHs 
clusterizing there. Section [VIII contains our conclusions. 



II. PBH MASS SPECTRUM CALCULATION 

A. General formula 

The calculation of PBH mass spectrum in Press- 
Schechter formalism [45j is based on the expressions 
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Here, the following notations are used: 5r is the initial 
density contrast smoothed on the comoving scale R, M is 
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the smoothing mass, a R {M) is the mean square deviation 
(the mass variance), 



oo 

°r(M) = Jv s (k)W 2 (kR) 



dk 

T' 



(3) 



Vs(k) is the power spectrum of primordial density per- 
turbations, W(kR) is the Fourier transform of the win- 
dow function (in this work we use the gaussian one, 
W(kR) — exp(— k 2 R 2 /2) ), pi is the initial energy den- 
sity. It is assumed that the process of reheating is very 
short in time, so, the end of inflation practically coincides 
with a start (at t = £j) of the radiation era. 

Fourier transform of the (comoving) density contrast 

is 



S k (t) 



k 
~a~H 



*k(«), 



(4) 



where is the Fourier transform of the Bardeen po- 
tential. In the approximation = const, we have, for 
the radiation dominated epoch, 8k ~ cl , that is the well- 
known result [49j . In this work, however, we will not use 
this approximation and, in opposite, we will explicitly 
take into account the time dependence of the Bardeen 
potential. 

The power spectrum of the density perturbations, cal- 
culated at some moment of time, is 



V s (k,t) 



-Xkrf 



(5) 



where r is the conformal time (r = (ai?) _1 for the radi- 
ation epoch). 

The comoving smoothing scale, R = l/k R , is con- 
nected with the smoothing mass M by the expression 



M \ - 2 / 3 fc? 
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where Mi, cii and Hi are the horizon mass, cosmic scale 
factor and Hubble parameter at the moment U. Eq. ([5]) 
immediately comes from two simple formulas: horizon 
mass at tj is given by 



47T , 

Mi = —tipu 
and the smoothing mass is 

4-7T af 

M = Tk% p - 

From these relations, one obtains 
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Mi 
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For the practical use of formulas ([T] [2]) one must ex- 
press the PBH mass Mqh through the smoothing mass 
M and the density contrast at the moment of the collapse 
(which is approximately equal to the density contrast at 
horizon crossing, 8 R ). 

It follows from Eq. ^ that the smoothed density con- 
trast at the initial moment of time, Sr, and at the time 
of horizon crossing, 5^ , are connected by the expression 



8 H 



$k-(th) 



(Th) 
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Here we use the approximation according to which the 
smoothed density contrast is proportional to its Fourier 
transform at some characteristic value of fc, fc = k*. If 
the power spectrum V$(k) increases monotonically with 
k, then, evidently, k* w ku : = 1/R. If the spectrum has 
a maximum at some value of fc, fc = kg, the reasonable 
estimate for fc* is 



fc* 



k R 
fco 



for k R < k Q , 
for fcfl > fco . 



(11) 



B. Gravitational collapse models 

The connection between values of the smoothing mass 
M, density contrast 5^ and PBH mass Mbh can be ex- 
pressed in the general form 



M BH = f(M,6^;Mi 



(12) 



The concrete expression for the function / depends on the 
model of the gravitational collapse. In the model of the 
standard spherically-symmetric collapse the connection 
is quite simple: 



M BH = (6%) 1 ' 2 M h 



(13) 



Here, Mh is the horizon mass at the moment of time, t = 
th, when regions of the comoving size R and smoothing 
mass M cross horizon. According to Carr and Hawking 
[50j. 1/3 < 8% < 1. The derivation of Eq. (JT3D is given 
in the Appendix [AJ 

The horizon mass at the moment when a region of the 
comoving size R crosses horizon is 



47T 

M h = —{aRfp(a) 



(14) 



(at this moment, a/kn = H : ). Using 0, one obtains 

(^H. t ) 2 



M h = MiHi— = Mi 
kR 



k R 



(15) 



Finally, using ©, the connection between Mh and M is 
derived 0: 



and putting Hi — l/ti we come to Eq. ©. 



M h = M} /3 M 2/3 . 



(16) 
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This equation connects horizon mass at the moment of 
time th when the perturbed region crosses horizon with 
the initial (at U) horizon mass and the initial (at ti) co- 
moving mass M. It is seen that the comoving mass de- 
creases with time (M/, < M) (see, e.g., |5l|). 

From (|13t I16|) one has the expression for the function 
/ for the Carr-Hawking collapse: 



This ratio is constant with time. It is clear from these 
arguments that a time dependence of ^ k (given, e.g., in 
Fig. [IJ will affect the form of the PBH mass spectrum 
(see Sec. EJ. 



f(M 5 ft' Mi) = (5 ft) 1 / 2 M 2 / 3 M 1 ^ 3 (17) Time dependence of the gravitational potential 



In the picture of the critical collapse [52| the corre- 
sponding function is 

f(M, 6% ; Mi) = k c (5i - 5 c )^M 2 / 3 Ml /3 , (18) 

where 5 C , j c and k c are model parameters. In this work 
we will accept the following set of parameters, which is 
in agreement with the recent calculations [53|, [54[ : 



0.45, 7 C = 0.36, k c 



(19) 



One should note that the neglect of a time dependence 
of the gravitational potential greatly simplifies a calcu- 
lation of the PBH mass spectrum because the smoothed 
density contrast 5r enters the PBH mass spectrum ex- 
pression only through the ratio 5r/<jr. In the limit of 
the constant ^ k one has, from Eqs. ([15]) and ©, 



For a calculation of a time dependence of the gravita- 
tional potential we use the approach suggested in works 
of Ref. [EH . We assume that reheating is rapid and, cor- 
respondingly, the moment of the end of inflation coincides 
with the moment of a beginning of the radiation epoch. 
The only difference from the consideration of j5f| is that 
we do not assume, in general, the validity of the slow-roll 
approximation at the period near the end of inflation. 

The general expression for the Fourier transform of the 
Bardeen potential at the radiation epoch is 



(rad) 



(t) = A k f A (x) + B k f B (x) , 



(23) 
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R —5 H 
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In the same limit, an expression for the mass variance is 
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It follows from Eqs. flU and (dU) that 
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f A (x) 



f B (x) 



a; cos x 



: (x sin x + cos a;) , 



(24) 



(25) 



where x = c s kr, c 2 s = p/p = 1/3 is the sound speed. The 
connection between ty k and Fourier components of the 
curvature perturbation 1Z in gauge invariant cosmological 
perturbation theory is given by two expressions (see, e.g., 
the review (56[): 



3#**(r) + — ^— **(r) = -(1 + w)K k (r) . (27) 

For obtaining the coefficients A k , B k we need two conditions for ^ k {r). The first one is the Eq. |26|) calculated for the 
end of inflation, when lZ k = TZ k 1 , ^ k = VP^ , t = n, and aH = aji2j. The second one is the Eq. (|2"T|) calculated 
for the beginning of the radiation era when lZ k = 7Z k ad \ ^ k = , and w — wq — 1/3. The junction conditions 

at the transition time are [[13], [II], [HH]: 



(28) 
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K 



( inf )( Ti ) = 74 rad) (r 4 ). 



(29) 



Using Eqs. l[26 | [27 | and 



one obtains for the coefficients A k , B k the following two conditions: 
A k f A {x t ) + B k f B ( Xl ) = (-^ 2 (1 + «,) ftj^Tj) , 



(30) 



[A k f A ( Xi ) + B h f B ( Xi )] 



3w n 



[A k f A ( Xi ) + B k f B {xi)] = -(1 + w )TZ k int) (n) 



(31) 



Two inputs in these expressions, 1Zj, ; (tj) and 72-^ (rj), can be determined by a numerical calculation (see, e.g., 
[3ll |) using the concrete inflation model. 

In the present paper, we assume, for simplicity, that the right-hand side (rhs) of Eq. (|30[) is equal to zero. It does 
not mean, however, that we exclude a consideration of non-slow-roll models because the rhs of Eq. pp[) contains, 
except of the term (1 + w), also the time derivative of TZ(k) which is typically small due to freezing out of TZ(k) outside 
of horizon (if the comoving size fc _1 crosses horizon some time before the end of inflation, see, e.g., [3l|). The general 
case when violations of slow-roll are continuing up to the end of inflation (and when, correspondingly, the rhs of Eq. 
([30| cannot be neglected) will be considered in a separate paper. 

As a consequence, ^^^(i") is proportional to 72-jjf (r^), and one has, instead of Eq. ([5]), the simple connection 
between density perturbation at any time and curvature perturbation at initial moment of time 



V s (k,t) 



where the expression for \E , / c (t) is given by [5f 



2TZ {iDi) (n) 

^k(i~) = — g — —[(x — Xi) cos(x — Xi) — (1 + xxi) sin(x — Xi 



(32) 



(33) 



Correspondingly, in this approximation only one in- 
put value, lZ k ini \ri), is needed for a calculation of the 
Bardeen potential. 

In Fig. Q] typical results of the calculation using Eqs. 
(|3"3")l and are shown. It is clearly seen from the figure 
that in our approximate approach (instantaneous transi- 
tion from inflation to the radiation era) the evolution of 
the Fourier transform of the density contrast starts from 
zero value. In practical calculations of PBH production 
we displace the starting moment from t% to t[ ignoring 
this production inside the short time interval At = tj — f,-. 
Specifically, the shift At was determined by the condition 
lg(f^/fj*) = 0.1, throughout all numerical calculations. 



D. Power spectrum with maximum 

One of scenarios in which significant PBH production 
is possible is the one with the power spectrum Vn(k) 
having a strong bump near some value k — kg in the 
region of small scales. Such models have been considered 
in many papers, including [H [H M, [H, HI, H HH . 

Having the exact form of the power spectrum Vn{k), it 
is possible to calculate the PBH mass distribution using 
the formalism introduced in the previous sections. How- 
ever, to be more general, it is convenient to introduce 



some parametrization of ■P-R.(fc), and use it to calculate 
PBH mass spectra. For example, we find that power 
spectra with peaks considered in [3 11 ] can be rather well 
parameterized in the region of the maximum with a dis- 
tribution of the form 



lgr n (k)=B+(lgV^-B)exp 



(lgfc/fco) 



2£ 2 



(34) 



(it is similar to the primordial spectrum parametrization 
used in [50l ] and in our previous work [60|). Here, the 
value of B is known from observations at large scales 
(B w —8.6), and its exact value does not affect the PBH 
production rate. V\ is the maximum value approached 
by the power spectrum, ko and E are parameters deter- 
mining the position and width of the peak. 

It is clear that any power spectrum with maximum can 
be more or less well approximated in the most important 
region (near the maximum) just by three parameters, de- 
termining the position of the peak, its height and width. 
Thus, the results from taking the spectrum in such a form 
will be quite general, because PBH production is mostly 
determined, of course, by the region of the spectrum near 
the maximum. 

The connection between the wave number k^ and hori- 
zon mass Mh at the moment when the mode enters the 
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FIG. 1: The dependencies of the gravitational potential 
(upper panel) and the density contrast 8k (lower panel) on 
time. H.C. denotes the moment of horizon entry for the mode. 
The wave number taken is k = Q.2t~ x , and normalization is 
7if f) (r i ) = 0.1. 
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lg M BH [g] 



FIG. 2: PBH mass spectra calculated for the models of the 
Carr-Hawking collapse (solid line) and the critical collapse 
(dashed line). The following set of the parameters was used: 
E = 3, = 0.02, Ml = 10 14 g. 
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16.0 16.5 17.0 17.5 18.0 

lg M BH [g] 

FIG. 3: PBH mass spectra for the standard collapse case 
with P£ = 0.0172, M£ = 10 17 g (for all curves); from bottom 
to top, E = 1, 3, 5. 

horizon is determined from ([6]) and (|16|) and is 

M h (k R ) = -2 , (35) 

k r 

and further we will use a notation MP = Mh(ko), which 
is the horizon mass corresponding to a value of fco. 

The result of PBH mass spectrum calculations with the 
particular set of parameters is shown in Fig. O It is seen 
that the mass spectrum strongly depends on the model 
of the gravitational collapse (at the same values of the 
fluctuation spectrum parameters). One should note that 
on this and some following figures we show on the vertical 
axis the quantity n BB X (ai/io) 3 j which is the comoving 
number density of PBHs, and is independent on the re- 
heating temperature Trh in the limit ki = CLiHi 3> fco- 
The case when values of fc ~ ki are important for the 
PBH production will be discussed in Sec. fVl 

Figure [3] shows a dependence of the PBH mass spec- 
trum on the parameter E in Eq. (13~3| . We can see that 
this dependence is rather strong, that comes from the 
fact that ubh has an exponential sensitivity to the mass 
variance an(M) (which, in turn, is determined by the 
shape of the spectrum, including its width). We also see 
that the mass spectrum is always rather strongly peaked 
near the PBH mass ~ Af°. 

The fraction of an energy density of the Universe con- 
tained in PBHs today, fipBH, assuming that a mass of 
the produced black hole does not change in time, is 

fipBH = — (—) I M BH n BH {M BH )dM BH (36) 
Pc\a J J 

(p c is the critical density). This formula is rather accu- 
rate for black holes with initial mass M BB M„ , where 
M* « (3i ao) 1/3 ~ 5 x 10 14 g is the initial mass of PBH 
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which reaches its final state of evaporation today [6l| . 
a Q = 8.42 x 10 25 g 3 s~ 1 , and t is the age of the Universe. 

One should note that for the concrete PBH mass spec- 
tra shown in Fig. [31 we have, for the case of E = 5, 
the inequality Qpbh > 1, which is an inconsistent re- 
sult, so the production of such a big amount of PBHs is 
forbidden. We will discuss the possibility of a constrain- 
ing the model parameters in the distribution of Eq. (134p 
in detail further, in Sec. El The value 7?£ = 0.0172 
chosen for the calculation of the spectra shown in Fig. 
|3J corresponds, in the case of E = 3, to the value VL PBH 
equal to « 0.23. This is in the range needed to explain 
the observed amount of non-baryonic dark matter [62| . 
n ilbm /i 2 = 0.106 ± 0.008 (with la uncertainty). For the 
case of E = 1, it appears that SIpbh ~ 10 -10 , probably 
too small of a fraction to produce observable cosmological 
consequences. 



III. NEUTRINO AND PHOTON SPECTRA 
FROM PBHS EVAPORATIONS 

PBH evaporation, predicted by Hawking [63[, causes 
evaporating black holes to produce an isotropic extra- 
galactic photon and neutrino backgrounds, which can be, 
at least in principle, measured experimentally. It can be 



calculated having the PBH mass spectrum, and the com- 
parison with observations can be made. 

Evolution of a PBH mass spectrum due to the evapo- 
ration leads to the approximate expression for this spec- 
trum at any moment of time: 

2 

n BH (m, t) = (3crf ™ m3)2/3 "fig (( 3ai + m3 ) V3 ) > 

where a accounts for the degrees of freedom of evaporated 
particles and, strictly speaking, is a function of a running 
value of the PBH mass m. In our numerical calculations 
we use the approximation 

a = const = a(M%£ x ), (38) 

where M^ x is the value of Mbh in the initial mass 
spectrum corresponding to a maximum of this spectrum. 
Special study shows that errors connected with such an 
approximation are rather small. In this work we use 
parametrization of the function cx(Mbh) presented in 

The expression for the extragalactic differential energy 
spectrum of neutrinos or photons (the total contribution 
of all black holes) integrated over time is [48| 



S(E) 



C I *7T (S)V nM af + m3 ) V3 ) •*W + *),m)e-*<*«>- m2dm 



4irJ a \a J J a " V ' J rv v " ' (3ai + m 3 ) 2 / 3 ~ 

(39) 

= j F(E,z)d\g(z + l). 



In this formula, (Zi, a, and ao are cosmic scale factors 
at ti, t and at present time, respectively, and ip(E,m) 
is a total instantaneous spectrum of the radiation (neu- 
trinos or photons) from an evaporation of the individual 
black hole. The exponential factor in Eq. (|39[) takes into 
account an absorption of the radiation during its prop- 
agation in space. The processes of neutrino absorption 
are considered, in a given context, in 48]. In the last line 
of (|3"9"]l we changed the variable t on red shift z using the 
flat cosmological model with J7a ^ for which 



i (ft M (i + z) 3 + tti?.(i + z) 4 + a0 



-1/2 



H 



(1 + *) 



(40) 



(H = 100/1 km/(s • Mpc), n R = ft 7 + n„), with the 
set of the basic cosmological parameters taken from [62j : 
h = 0.73, Q M h 2 = n nbm h 2 + n B h 2 = 0.128, Oy/i 2 = 
2.47 x 10~ 5 , Q\ = 1 — Q M — Ct R . The neutrino fraction 
in the total energy density, Q^, is calculated assuming 
that the neutrino masses are negligible, in this case = 



I(ti) 4 ^ 3 x 3 ^7 ~ 0-68S! 7 (the term with fi^ is essential 
in Eq. gDJ) only for z > Qm/^r)- 




FIG. 4: Absorption cross section a a for photons, in units of 
7rr^. Dashed line is the result of the numerical calculation 



[66| , solid line is the approximation (|4"3| . 
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The direct Hawking flux of evaporated particles (per 
degree of freedom) is 



<p(E,M BH ) 



dN 



dEdt 2tt exp(E/T H ) - (-1) 2 



(41) 



where T s is a grey-body factor (T s = a s E 2 /tt for massless 
particles, where a s is the absorption cross section), and 
Tjj is the Hawking temperature 



Ti 



H 



1.06 ff^ GeV. 
%-kMbh \MbhJ 



(42) 



In the high energy limit (E ^> Th), the cross section a s 
approaches a constant value independent of s, and T s in 
this case is 



27 E 2 
64^ if' 



(43) 



For lower energies, however, this factor needs to be cal- 
culated numerically. 

For photons, we use the following approximate 
parametrization of this function 64 1 



(exp[9.08 - 1.71(£ 7 /T H )] + 1) 



4 ( E -< V 
L 3 \^47rT H ) 



E-t > 2.5T H , (44) 



E 1 < 2.5T H 



Here, the low-energy part comes from analytically ob- 
tained relation a s =i(E — * 0) = ^nr 2 (Er g ) 2 [6l[ (r g is 
the Schwarzschild radius). The higher energy part is just 
an approximation smoothing the oscillations in a s , see 



Fig. [U The form of Eq. (|44|) at high energies is analo- 
gous to one used in [65j for neutrinos, 



r„ = 



27 E 2 
64^2 if 



0.075 



0.925 



exp[5 - 1.607(E/T H )} + 1 



(45) 

The instantaneous spectra of neutrinos and photons 
from PBH evaporations were calculated using the model 



of Ref. [66j and parametrization of Ref. 67]. In these 
papers the particles directly emitted from the black hole 
(through the Hawking mechanism) as well as those pro- 
duced in processes of quark fragmentations and subse- 
quent decays of mesons had been taken into account. 

In Fig. [5] the red shift distributions of the differential 
energy spectra S(E) are shown for the neutrino or photon 
with energy 1 GeV. One can see that the distributions 
are more wide in the case of the critical collapse (because 
the PBH mass spectrum in this case has a long "tail" of 
small masses). The absorption factor e~ r is efficient: 
in case of photons, the maximum value of z for which 
the absorption can be neglected is around 700 (68|; for 
neutrinos the absorption is important at z > 10 6 -j- 10 7 

Hi. 



> 

o 
O 
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bo 
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FIG. 5: Red shift distribution of the integrand F(E, z). Thick 
lines are for neutrinos, thin lines are for photons (absorption of 
7-rays at z > 700 is not shown in this figure). Dashed curves 
represent the case of the critical collapse, solid ones corre- 
spond to the standard collapse. PBH mass spectra shown in 
Fig. [2] were used in the calculation, E — lGeV for all cases. 



IV. CONSTRAINTS ON THE POWER 
SPECTRUM WITH MAXIMUM 

Assuming that Vfi(k) has the form of Eq. (f3~4l -. we can 
calculate the produced PBH mass spectra, compare the 
possible consequences of the existence of such amount of 
PBHs with observations, and put limits on the param- 
eters of this distribution. The calculation of the actual 
constraints in this paper (as well as in our previous work 
[rJol ]) uses the following basic observational facts. 

1. The differential energy spectrum of the extragalac- 
tic photon background in the wide region of energies is 
known (69| . For our purpose, the most interesting range 
is about E 1 ~ 1 MeV GeV, and the order of mag- 
nitude of the flux is roughly given by E 2 S 1 (E) ~ 10~ 6 
GeV cm -2 s _1 sr _1 in this region. We obtain our con- 
straints just from the condition that the diffuse gamma 
ray flux produced by PBHs does not exceed the observed 
extragalactic one. 

2. According to the data of Super-Kamiokande exper- 
iment (7fjj . the electron antineutrino background flux in 
extragalactic space is constrained by the inequality 



<f>(E Pe > 19.3 MeV) < 1.2 cm 



-v 1 . 



(46) 



We integrate the calculated flux and put constraints from 
the condition that it does not exceed this limit. We note, 
however, that the limit (|46|) is not a model-independent 
one. It was obtained assuming exponentially falling neu- 
trino spectrum, which is an expectation for the extra- 
galactic diffuse flux from supernovae. The neutrino flux 
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FIG. 6: Constraints on the parameter (giving the height 
of the peak in Eq. (I34|l '). that were obtained in this work 
(for this figure, we used E = 3). The forbidden values of 
the parameter V%_ lie above the curves shown in the figure. 
Constraints following from the neutrino and the gamma ray 
experiments are shown together with the energy density con- 
straint. Dashed lines correspond to the model of the critical 
collapse, solid lines represent the results obtained using the 
standard collapse picture. 



from PBHs in the mass region we consider is falling as 
E~ 3 , which, in fact, will change the constraint (|4"(51) . Es- 
timations show, however, that this change is about a fac- 
tor of 2, and it is not very significant for putting con- 
straints on such parameters as V^. Similar uncertainty 
also comes from the inclusion of neutrino oscillations in 
the analysis, which we do not take into account. 

3. The fraction of the energy density of the Universe 
contained in PBHs, f2pBH, which can be calculated us- 
ing formula (|36[) . cannot exceed the one for non-baryonic 
dark matter, fi n bm [62| . This constraint is important for 
black holes with initial mass Mbh > M*, i-e., ones that 
did not evaporate up to the present time. 

The actual condition we use is flpsH < 0.3 (having 
all the uncertainties considered, f2 n bm = 0.3 is still al- 
lowed by the observations within the 3ct confidence in- 
terval, however, the uncertainty in this number does not 
seriously affect the power spectrum constraint). 

The resulting constraints obtained from all the condi- 
tions discussed are shown in Fig. [75] The value of the pa- 
rameter E characterizing the width of the gaussian distri- 
bution in Eq. (f3~4"|) was fixed throughout all calculations 
(E = 3). The constraints are given as a function of the 
horizon mass MP. The connection of M® with the corre- 
sponding wave number fco is very simple and is given by 
the formula following from Eq. 



3 x 10 23 , r 
fc w Mpc~ 



(47) 



It follows from these results that the constraints are 
stronger under assumption of the standard Carr-Hawking 
collapse, which is expected because in this case the 
threshold density contrast leading to PBH formation is 
smaller. The constraints are independent on the reheat- 
ing temperature if it is high enough to provide the condi- 
tion Mi < M£ (this is true for all PBH masses considered 
if T RH > 10 11 GcV). 

The second important result is that the constraints 
based on the neutrino emission of PBHs are compara- 
ble with those following from the photon emission and 
from the gravitational constraint. At the region of small 
horizon masses (large fco), where large red shifts are im- 
portant, the constraints from the neutrino emission are 
stronger. 

Fig. El is analogous to one presented in our previous 
work [60|| . In the present paper we add the gravitational 
constraint on Qpbh and more carefully investigate the 
mass region M® ~ fO 16 -j- 10 17 g. In this region, con- 
straints from energy density and from gamma ray back- 
ground coexist. We do not extend in the present paper 
the neutrino constraints to masses M% > M* because the 
neutrino background in this case strongly deviates from 
the simple behavior ~ E~ 3 and additional analysis is 
needed to estimate how the limit of Eq. ([46]) will change 
in such a situation. 



CONSTRAINTS ON THE RUNNING MASS 
MODEL FROM PBHS 

The running mass inflation model was proposed in 
I and further studied in many papers including 
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Z3, IZJ, lZ5|, IZfil, IZ7J- The model predicts a rather strong 
scale dependence of the spectral index, possibly allowing 
large values of Vtz(M1 a ^ small scales, which can lead to 
PBH production [[78[, [zD], [Hj]]. Th e first constraints 
on the running mass model parameters from PBHs were 
obtained in works [78|, [7t| . In previous works the PBH 
mass spectrum in the running mass model have not been 
calculated in detail. We do it in the present paper to 
demonstrate how the general formalism developed in Sec. 
HI1 works. 

We will also show how it is possible to constrain such 
observable parameters as the spectral index measured at 
cosmological scales, no, and its running n' , assuming the 
running mass model is correct. 

The potential of the running mass model takes into 
account quantum corrections in the context of the softly 
broken global supersymmetry and is given by the formula 



V = V a 



1 



to 2 (In 



(48) 
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The dependence of the inflaton mass m(ln <f>) on renor- 
malization scale <f> is determined by the renormalization 
group equation (RGE). 

Apart from Vq, the shape of the running mass potential 
is determined by parameters c and s, which are connected 
to the other quantities by 



Vo 

'A/2 



dm 2 



din i 
Hi 



4>=<po 



(49) 
(50) 



(H'j = Vo/3Mp, and 0o is the value of the inflaton field 
corresponding to the cosmological scale ko). They are 
also linked to observables no and n' : 



n 



1 « 2(s — c) , n' « 2sc , 



(51) 



and we see that measuring uq and n' allows us to recon- 
struct the potential shape in this model, and the only 
free parameter left is Vo (or Hi). From the theoretical 
point of view [77| , Hi can lie in the wide range of values 
from, say, Hi ~ 10 4 GeV for anomaly-mediation case to 
Hi ~ 10 -3 GeV for gauge-mediation. Assuming instant 
reheating at some value of the inflaton field, from the 
energy conservation condition we have: 



30 9 * 1rh 



(52) 



(<7* ~ 100), and the reheat temperature turns out to be 
in the region 10 7 -j- 10 11 GeV for the given values of Hi. 
For further calculation we will use two values of Tr#, 
namely, 10 8 and 10 10 GeV. 

It was noted in [3l|, that for the running mass model, 
the numerical calculation of the power spectrum is nec- 
essary in the region of large k values. In this work we 




lgk[Mpc _1 ] 



FIG. 7: Power spectrum Vn{k), calculated for the running 
mass model, with Trh = 10 10 GeV, no = 0.97, from top to 
bottom, 10 3 n = 4.3,4.5,4.7 . 
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FIG. 8: PBH mass spectra in the running mass model, cal- 
culated for the following set of parameters: Trh = 10 8 GeV, 
n = 0.96, n' = 5.33 x 10~ 3 (upper panel); T RH = 10 10 GeV, 
no = 0.96, n' = 4.4 x 10~ 3 (lower panel). Solid line rep- 
resents the result for the standard collapse case, dashed line 
corresponds to the critical collapse. 



also use the numerical calculation of power spectra, not 
relying on the analytical approximations. Several results 
for Vn(k) are given in Fig. [3 which shows how sensitive 
is the value of Vn near the end of inflation to the change 
in n' . 

Some examples of the PBH mass spectrum calculation 
are given in Fig. [8] for two models of the collapse con- 
sidered above. The order of magnitude of the PBH mass 
produced is equal to the initial horizon mass Mj, which 
is given by 



M, 



0.038- 



m 



Pi 



1/2 



T 2 



(53) 



The real PBH mass spectrum, and, particularly, its max- 
imum, is, however, determined by the dynamics of their 
formation, including the dependence of the gravitational 
potential on time. This dependence is very important 
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FIG. 9: Constraints on the running mass model coming from 
PBHs, for two values of parameter Trh- The forbidden re- 
gions of parameters lie above the corresponding curves. Solid 
thick lines are constraints for the standard collapse model, 
dashed thick lines are those for the critical collapse. Thin 
lines show, for a comparison, the values of no and n' , for 
which Vn = W 1 (dotted line) and Tn = 10~ 2 (dashed line) 
at maximum value of k, k = ki- 



rn the region of smallest possible PBH masses produced, 
i.e., in the case Ur ~ ki. For example, in the model 
of the standard collapse with a monotonous power spec- 
trum we would expect most of PBHs to be produced near 
the smallest possible mass, i.e., near Mrh ~ ^/l/3 M, 
(e.g., M BH ~ 7 x 10 14 g for T RH = 10 8 GeV). From 
Fig. [8] we see, however, that the actual maximum is near 
Mbh ~ 3 x 10 15 g in this case, which is due to suppres- 
sion of PBH production at smaller masses. 

Fig. [HI shows the constraints on (no,n' ) obtained as- 
suming that the running mass model is correct. Recon- 
structing the potential from values of no, n' and Trh , we 
calculate the power spectrum Vfi(k) in the whole region 
of wave numbers up to ki, and the corresponding PBH 
mass spectrum. The method of obtaining the constraints 
depends on the characteristic PBH mass produced: for 
Trh = 10 10 GeV, the PBH mass in maximum of the dis- 
tribution is about 3 x 10 11 g (see Fig. [5]), and neutrino 
background data was used to constrain such models, be- 
cause, as we have seen in Sec. IIV1 neutrino experiments 
arc adequate to constrain PBHs of such initial masses. 
For Trh = 10 8 GeV, the PBH mass in maximum of the 
spectrum is about 3x 10 15 g, and constraints are obtained 
using the extragalactic gamma ray background. 

We see from Fig. [3]that, assuming running mass model 
is really responsible for inflation and production of per- 
turbation spectrum, running of the spectral index at cos- 
mological scales is strongly constrained by PBH overpro- 
duction and cannot be larger than (3-^6) x 10 -3 (the 
exact value of the constraint depends on other parame- 



VI. DARK MATTER AND PBHS 



PBHs produced in the early Universe, if they are heavy 
enough not to evaporate till the present day, can con- 
tribute to dark matter, and potentially it is possible to 
explain all the non-baryonic dark matter with them. For 
example, the mass spectrum shown in Fig. [3] (for the 
case of £ = 3) give s Qpbh ~ 0.23. 

In works [43, [44J , PBHs clustered in the Galactic bulge 
are proposed to explain the observed 511 keV photon line 
from the Galactic center [§3, HH . Positrons evaporated 
from PBHs annihilate with electrons from the interstellar 
medium thus producing the observed photons. Assuming 
that Solar System lies at a distance of ro = 8.5 kpc from 
the Galactic center, the total positron production rate 
needs to be about 3 x 10 43 s _1 , and the corresponding 
characteristic value of the PBH mass should be in the 
range ~ 10 16 -j- 10 17 g (this was shown in ref. [44| from 
calculations of the photon background in extragalactic 
space and the photon flux from the Galactic center from 
PBH evaporations). We will show in this paper that 
our distribution presented in Fig. [3] (for £ = 3) is also 
capable to explain the needed positron production rate. 

Assuming that the PBH mass distribution in Galaxy 
does not change compared to the initial one, for the PBH 
mass distribution in the Galactic center Nbh we can 
write a simple proportionality relation 

Nbh(Mbh) = C n BH {M BH ), (54) 

where C will be determined by the total mass of PBHs: 

J N BH {M BH )M B HdM BH = Mtot- (55) 
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FIG. 10: Gamma ray flux from PBHs clustered in the Galactic 
center, assuming their total mass is 1O 9 M0 and mass distri- 
bution is the same as in Fig[3j(E = 3) (thick line). The thin 
line shows the measured background flux. 
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The gamma ray flux on the Earth surface produced by 
such PBHs is 



dN., 



dE^dSdt 



(56) 



N B h{M BH ) 



dN 1 
dE-ydt 



(M BH ,E 7 )dM BH 



and the positron production rate in the Galactic center 
is 



dN e+ 
dt 



(57) 



J J N BH (M BH ) (M BH , E e+ )dM BH dE e 



The total mass needed can be expressed through the 
observable positron production rate and the initial mass 
spectrum by the formula 



dN 



dt 



J n BH (M BH )M BH dM BH 



J n BH (M BH ) dE'+dt ( M bh ,E e +) dE e + dM BH 



(58) 



In our case, we obtain M tot ~ 10 9 M Q . If these black 
holes are uniformly distributed in the spherical region of 
radius r ~ 600 pc, the clusterization factor needed for 
that is about 



c 



M t , 



|7rr 3 Vt PBH p c 



3 x 10' 



(59) 



The gamma ray flux from such PBHs calculated us- 
ing Eq. (j56"|) is presented in Fig. [TU1 together with the 
measured continuum, given by the approximate formula 



d^ r 



dE^ 



7 GeV -i cm-V 



E-. 



511 keV 



-1.75 



(60) 



and we see that the flux from PBHs does not exceed the 
observational one. 

For the calculation of diffuse extragalactic photon flux, 
clusterization of PBHs is not important, because this flux 
is sensitive to the average number of PBHs and their 
mass distribution only. Thus, formula (|39[) is suitable 
for the calculation. It can, however, be much simplified 
because such PBHs (with mass sufficiently larger than 
M*) do not lose the significant part of their mass in course 
of the evaporation, and approximation M BB {t) ~ const 
is adequate, so we can neglect the quantity Sat in it. 
The absorption factor can also be dropped because the 
effectively working red shifts in this task are small. The 
diffuse gamma ray flux from our PBH mass distribution 
(Fig. [SJ £ = 3) is shown in Fig. [TT] We see that, again, 
the experimentally measured flux, given by 0, H^] 



d$ c 



6.4 GeV^cm-W- 1 



1 MeV 



-2.38 



is not exceeded. 
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FIG. 11: Diffuse gamma ray flux produced by PBHs with the 
mass spectrum shown in Fig [3] (E = 3) (thick line) and the 
measured extragalactic background (thin line). 



VII. CONCLUSIONS 

1. We have developed a method of PBH mass spec- 
trum calculation, which is based on the Press-Schechter 
formalism and takes into account the dependence of the 
gravitational potential on time. We found that it is more 
convenient and natural to use in the formula for the 
mass variance (Eq. ([3])) the gravitational potential $&(£) 
calculated as a function of time rather than the trans- 
fer function T(k,i) (the latter is defined, e.g., in [4l|). 
To derive the PBH mass spectrum, the power spectrum 
of curvature perturbations calculated (using inflationary 
models) at the end of inflation is needed as an input. 

2. PBH mass spectra are calculated in two inflation- 
ary models predicting large amplitudes of the curvature 
perturbation spectrum at some wave number k, namely, 
in the model with the power spectrum with a peak and in 
the running mass model. Examples of such a calculation 
are given in Figs. [21 El and [HI 

3. Constraints on parameters of both these models are 
obtained using available data on the extragalactic pho- 
ton diffuse background and the upper limit on the ex- 
tragalactic neutrino background, obtained in the Super- 
Kamiokande neutrino experiment. 

The results are shown in Figs. [6] and 03 It is seen 
that for the power spectrum with a peak, the height and 
the width of the peak can be constrained in rather wide 
region of wave numbers determining its position in the 
spectrum. For the running mass model we obtain the 
constraint on the spectral index running at cosmological 
scales, n' < (3 -j- 6) x 10 -3 (the exact constraint value 
depends on no and T BB )- 

4. It follows from Fig. [S] that constraints obtained 
from neutrino experiments are stronger than those ob- 
tained from extragalactic photon background, in the re- 

(61) gion of horizon masses from ~ 10 11 to ~ 10 13 g. So, 
constraints from neutrino experiments and from measure- 
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ments of diffuse photon background are complementary. 

5. PBHs can be the main contributor to the dark mat- 
ter in the Universe. It is shown that if, in this case, the 
PBH mass distribution is peaked near ~8x 10 16 g, as in 
Fig. [3] with £ = 3, the experimentally observed galactic 
511 keV photon line can, in fact, be caused by annihilat- 
ing positrons produced by evaporating PBHs clustered 
in the Galactic center. The total PBH mass needed for 
this to happen must be about 10 9 M Q , and clusterization 
factor needed is C ~ 3 x 10 7 . 
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APPENDIX A: FORMATION OF PBHS IN A 
RADIATION-DOMINATED UNIVERSE 

In this Appendix we give basic formulas of the PBH 
formation process following closely the work |l8|. 

Evolution of a spherically symmetric region with den- 
sity p greater than the background density p is governed 
by the positive curvature Friedmann equation, 



H 2 (t) 



3m 



8?r 

—P(t) 



pi 



R 2 (t)' 



(Al) 



Here, H, p and R are Hubble parameter, density and a 
physical size of the region, respectively. In a background 
space the evolution equation is, correspondingly, 



H 2 (t) = 



8tt 
3m 2 pl 



p(t). 



(A2) 



An expansion of the overdensed region goes on accord- 
ing to the law [82] 



R~ t 



1/2 



(A3) 



and stops at time t c . This time is determined from the 
condition H(t c ) = and is given by the expressions 



t c . = S- U 



Si = 



Pi - Pi 



(A4) 



where U is the initial moment of time (which is chosen in 
such a way that H{ti) = H{ti) = Hi [18|). This ti should 
not be confused with t t of Sec. [TTJ The size of the region 
at t c is 



R(i c ) = R c = S\ 1/2 R(t t ) = 5 i i/z R,. (A5) 
The perturbed region will contract and, eventually, col- 
lapse, if it contains enough matter to get over pressure 
forces, i.e., if its radius exceeds the Jeans length, 

Rc ^ Rjcans — C s t c — j=t c . (-^6) 



-1/2, 



The corresponding requirement for the PBH formation is 



(A7) 



^/^C r^j Rc rS t c 



(the value of the upper bound is explained in [[50(, [37 
Dividing (|A7|) by t c , we obtain the expression for R c /i c , 



Rc _ Rj ^1/2 

tc t{ 



(A8) 



which scales, approximately, like a constant (independent 
on time) because 



■Rj ^° ^ft~i , Si 00 t{ 



(A9) 



Evaluating (|A8|) at horizon crossing (where R# = tn), 
we obtain from (|A7|) the necessary condition for the col- 
lapse (and the PBH formation) 



<5g<l. 



(A10) 



Here we used the notation for the density contrast at 
horizon crossing introduced in Sec. ITT1 

The mass of a black hole formed in a process of the 
collapse is given by 



47T 

M BH = yRcP(tc), 



where, from JA8D and (TA4 



R c = ^! /2 ic ; i c = W 1 ; 

Li 



P{tc) = Pt 



(All) 



(A12) 



(A13) 



Substituting (jAl~2|) and (fAT3|) to ([ATT]) , one obtains 



4?r 

M BH = -^r 



Ri r 

— 6 

t. 



1/2 



(A14) 



The whole expression (IA14|) . in its r.h.s., scales, ap- 
proximately, like a constant with time, similar with r.h.s. 
of Eq. (|A8p . Evaluating it at horizon crossing (R# = 
tn), we obtain the connection between black hole mass, 
horizon mass Mh and 5^ , given in Sec. |IIJ 



4-7T . 



Mbh - ^{^) 1/2 t H p H = M h {5^fl 2 . 



(A15) 



Using the connection between Mh , M and Mi given in 
Eq. (fi"6| of Sec. [Til one obtains the well-known formula 
for M BH (see, e.g., \M,M M)- 

M BH = (Si flHl^M^, S% « i = 7 . (A16) 
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